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Abstract. In this article, I briefly review the status of infrared effects which occur 
when using inflationary models to calculate initial conditions for a subsequent hot, 
dense plasma phase. Three types of divergence have been identified in the literature: 
secular, "time-dependent" logarithms, which grow with time spent outside the horizon; 
"box-cutoff" logarithms, which encode a dependence on the infrared cutoff when 
>■• ' calculating in a finite-size box; and "quantum" logarithms, which depend on the 

<~i , ratio of a scale characterizing new physics to the scale of whatever process is under 

Mh' consideration, and whose interpretation is the same as conventional field theory. I 

Q . review the calculations in which these divergences appear, and discuss the methods 

.y ' which have been developed to deal with them. 
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1. Introduction 

Looking out at the universe from the vantage point of Earth, we see a small fluctuation 
ST in the temperature of the cosmic microwave background radiation (CMB). Over 
the last decade, advances in observational astronomy have allowed us to commence a 
detailed study of (5T's statistical properties. We aim to compare these properties with 
those of small fluctuations which, according to our present ideas, are expected to have 
existed in the very early universe. 

These fluctuations evolve in a hot, dense plasma of tightly coupled baryons and 
photons. If we suppose the plasma era was preceded by an epoch of primordial inflation, 
then initial conditions for these fluctuations can be calculated from the parameters of 
the inflationary model. The variance predicted by this approach was computed in a 
now-classic series of early papers [H El El HI El E] , and later extended to the skewness 
m El El [ini im [12] and kurtosis [13 [12 IISI [16] . With sufficiently precise experimental 
data in hand, we can work backwards from observations of these statistical properties to 
determine what the initial conditions must have been. The result is a realistic prospect 
of constraining parameters in certain inflationary models. 

This ambitious programme was begun in earnest following the arrival of high-quality 
CMB maps made by the Wilkinson Microwave Anisotropy Probe (WMAP) [17J. Even 
more precise data are expected from the Planck Surveyor satellite. Once it had become 
clear that the statistical properties of 6T carried important information about whatever 
physics was operative during inflation, theorists were soon tempted to refine their 
calculations in the hope that additional detail could be extracted. A similar programme 
in electroweak physics had provided important clues about the unknown details of 
electroweak symmetry breakingqj and it was reasonable to explore the possibility of 
similar gains in cosmology. The observable statistical properties of the infiationary 
density perturbation are encoded in its ra-point correlation functions, with the leading 
contribution to each function typically suppressed by (n— 1) powers of the small quantity 
(iJ/Mp)^ [19], where H is the Hubble rate during infiation and Mp = (SvrG)"^/^ is the 
Planck mass. Two types of refinement were possible: either to calculate to lowest order 
in [H/M-pY at progressively larger n, or to calculate higher order corrections with n 
fixed. 

It had been known for a long time that perturbation theory beyond leading order 
in [H/M-pY was complicated by troublesome infrared behaviour. Sasaki, Suzuki, 
Yamamoto & Yokoyama considered a A0" theory with n ^ 3, coupled non-minimally to 
gravity, and evaluated each n-point function perturbatively [20l[21]. Sasaki et al. noticed 
that the n-point functions evaluated in de Sitter space had a potential divergence in the 
far future, which they called "superexpansionary." We will rephrase the analysis of 
Sasaki et al. in modern notation and discuss the meaning of these divergences in §2.11 
and P 

f More details of this "precision electroweak programme" can be found in the review Electroweak model 
and constraints on new pfiysics, prepared by the Particle Data Group [18] . 
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Later, Mukhanov, Abramo & Brandenberger studied a different sort of infrared 
effect — the accumulation of long wavelength fluctuations in an expanding universe 
[221 [23] . (See also Abramo & Woodard, Ref. [21].) This differed from the analysis 
of Sasaki et al, which had been based on a purely geometrical effect. Whether such 
divergences are physical is complicated by the question of gauge, initially studied by 
Unruh [25]. Later, Losic & Unruh were able to give a gauge- invariant argument that 
infrared terms generically become large in spacetimes close to de Sitter [261 l2ll ISHl 129] . 
We will discuss corrections of this type in §2.21 and §H Divergences of a third kind 
had long been studied by Prokopec, Tsamis, Woodard and their collaborators. These 
authors took the quantum nature of correlation functions seriously, and calculated loop 
corrections just as one would when studying scattering experiments. Unfortunately, 
the complexity of both the rules for calculating these loops and the computations 
themselves are somewhat greater than in Minkowski space. These calculations have 
been reported in a large literature. Together with quantum effects studied by 
others, including Boyanovsky, de Vega & Sanchez, a recent sample can be found in 
Refs. [30l[3ll|33[33l|3l[35l[36ll37l[38ll39lllQlliT]. We will consider effects of this type 
in 3231 

Perturbative calculations of correlation functions from inflation had been studied in 
the 1980s by Allen, Grinstein & Wise [7] and later by Falk, Rangarajan & Srednicki [8]. 
The sophistication of such calculations increased dramatically in the years following 
Maldacena's calculation of the three-point function in single-field slow-roll inflation 
[To] . At least in part, this increase was a consequence of the more complicated models 
of inflation which had been developed during the 1990s and early 2000s, and which 
theorists hoped would yield a distinctive pattern of correlations [42l [121 HHl [IB] ■ These 
developments are reviewed pedagogically in a recent article by Chen ^^. 

Since 2002, significant theoretical effort has been expended in refining our 
understanding of both inflationary correlation functions themselves, and the models 
which can generate significantly non-gaussian statistics. The calculations reviewed 
in this paper were a consequence of the reinvigoration of interest in infrared issues 
which followed this effort. Logarithms exhibiting secular growth with time were 
discovered in calculations of n-point correlation functions with n ^ 3 [SI |45l |46l I47ll48] . 
Other calculations, if taken at face value, gave divergent answers in the infinite- 
volume limit. Divergences could be avoided by carrying out the calculation in a finite 
box, but any such procedure left behind residual logarithms involving the box cutoff 
^9\ [sot B6l W7\ [51] |52| 153] . Further logarithms were encountered while studying the 
influence of unknown physics at energy scales even larger than the inflationary scale 
[5a[55l[56l[571|58l[59l|60l[6T]. 

Infrared effects are not a phenomenon unique to de Sitter calculations, but are 
present in many applications of quantum field theory. Their appearance typically signals 
the presence of a non-trivial background, which cannot be described by asymptotic 
in- and out-states containing a definite number of particles. In scattering calculations, 
infrared divergences allow putative fixed-particle states to develop a slowly varying field. 
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made out of an accumulation of arbitrarily many soft particles radiated on approach to 
or recoil from the scattering event. This "Bloch-Nordsieck" or "initial state radiation" 
phenomenon is studied in many textbooks on quantum field theory [62]. A similar 
effect can occur in the confined phase of quantum chromo dynamics, where soft particles 
can be radiated by partons moving inside nuclei. These particles are trapped by the 
strong QCD force, and change the background colour field. This effect may dramatically 
alter the mix of partons observed by an impinging probe such as a photon, and cannot 
be neglected in an accurate comparison with experiment |63l [Ml ESI IMl ETJ EH] . The 
effective W approximation is a related example, in which the same infrared effects 
control the W^ content of colliding nuclei. The results are very practical. For example, 
this method can be used to predict the rate of Higgs production from W^W~ fusion 

Throughout this article, Planck's constant h and the speed of light c are set to unity. 
Many expressions are written in terms of the reduced Planck mass, Mp ^ = SvrG. For 
brevity, a scalar field which contributes to the energy-momentum tensor Tab is referred 
to as 'active,' whereas a field making no contribution to Tab is referred to as a 'spectator.' 
In a cosmological context, such scalars are isocurvature modes. 

2. A zoology of logarithms 

In this section, I briefly review the sources of logarithmic divergences in inflationary 
correlation functions. 

2.1. Time- dependent logarithms 

A decade after calculations of the 2-point function had determined the variance in 
5T expected from an inflationary initial condition, the corresponding skewness was 
calculated by Falk, Rangaraian & Srednicki [8] in an approximation where gravitational 
interactions were neglectedljj At that time, data from the COBE satellite [71] had 
established the existence of temperature fluctuations in the microwave background at 
the level 5T ~ lO^^T, where T k. 2.75 K was the average microwave temperature over 
the sky. In inflationary models 5T /T is closely related to the parameter [H/M-pY . Since 
the skewness is proportional to [H/Mp)'^, it was already evident that any non-gaussianity 
would be very small. 

Let us suppose inflation is driven by a single scalar field with potential V. We 
define the slow-roll parameters by 

f It was later argued by Gangui et al. [S] and in more detail by Maldacena [TU] that gravitational 
interactions in fact provide a dominant contribution to the skewness, but the contribution calculated 
by Falk et al. is still present, if subdominant. See Ref. [44] or the review of inflationary non-gaussianities 
by Koyama in this volume |70| . 
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and use a conformal time variable, r, related to cosmic time t by r = J dt'/a(t'). 
When measured in this variable, horizon crossing occurs for a mode of wavelength k 
when kr = —1, and after N e-folds outside the horizon we have kr = —e~^. The 
infinite future of de Sitter space corresponds to |A;r| — )• 0. In this notation the result of 
Falk, Rangarajan & Srednicki can be written 

(<50(ki)50(k2)50(k3)). = {2nY6{k, + k2 + ka)^^^^ In \kn\ + ■ ■ ■ , (2) 

where a subscript '*' denotes evaluation at the time r^:, and '■ ■ ■' denotes terms of 
lower order in the slow-roll expansion, together with other terms which do not grow 
as I kr^ I — )■ [ISl HH] g| Higher powers of In | fcr* | may be generated at higher order 
in the slow-roll expansion, either by retaining such terms explicitly in the tree-level 
calculation or by including loop corrections, as we shall shortly explain. A similar effect 
can be observed in the massless limit of the calculation by Chen & Wang [73] . 

The result is a power series in In | A;r* | with an apparent divergence in the limit 
r* — )■ 0, as Sasaki et al. had predicted [20] . Series representations of this type 
were later applied to inflationary correlation functions by Gong & Stewart [711 ESj, 
who obtained them by systematically solving Mukhanov's equation [5] using a Green's 
function approach. Gong & Stewart remarked that the correction linear in logarithms 
was typically proportional to ~ e In |/i;r*|, whereas the correction quadratic in logarithms 
was typically proportional to ~ (eln |/i;r*|)^, and so on|§| One could therefore expect the 
series to break down as a predictive instrument when | In \kT^\\ ~ e^^. This occurs when 
the mode k is of order e~^ e-folds outside the horizon. 

Taking into account the construction of Gong & Stewart, and bearing in mind 
that T^ is the time at which we wish to evaluate the correlation functions, it seems 
clear that there is nothing mysterious about the appearance of powers of In |/cr*|. Such 
terms merely express that correlations evolve outside the horizon. This evolution is 
mainly driven by the evolving classical background field [45], but we shall see in §3.31 
that intrinsically quantum contributions can also be present. We conclude that, as in 
other applications of quantum field theory, properties of the background play a role in 

I If our purpose were only to calculate the density fluctuation produced by a single-field model of 
inflation, then the logarithmic divergence in Eq. ([2]) would be entirely spurious [25j . The density 
fluctuation is determined by the comoving curvature perturbation, C, which is conserved on superhorizon 
scales, C = 0, in the absence of entropy modes. Occasionally it has been argued that an estimate of the 
non-gaussian yield at the end of inflation can be obtained by setting Inj/cr*! ~ 60, but this example 
shows clearly that unless care is taken to account for possible cancellations any such procedure can be 
quite misleading. The possibility of such cancellations when reasoning with cutoff-dependent terms in 
a general effective field theory was pointed out by Burgess & London J72J. In the inflationary context, 
the pitfalls of this procedure were noted in Ref . [48] . I would like to thank Xingang Chen for discussions 
on this point. 

§ There is no general proof that the power series in ln|A;r*| always orders itself in this form, and 
according to a general theorem due to Weinberg (to be discussed in ^3.1\i there may be cases where it 
does not. In many practical examples where the slow-roll approximation applies, however, this structure 
appears. In these expressions, the symbol '~' is used to mean that the coefficient of In |fcr*| is a quantity 
of order e (but not necessarily e itself), wheras the coefficient of (In |A;t*|)^ is of order e^, and so on. 
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interpreting the physics of infrared divergences. In the present case, correlation functions 
can be evaluated at any time of interest provided the series of logarithms can be summed 
or ignored. Therefore, for observational purposes the key question is the value of e when 
modes of interest leave the horizon. We would typically wish to evaluate correlation 
functions at the end of inflation, where initial conditions for the subsequent evolution 
must be set. If e ^ 10~^, the end of inflation may be uncomfortably close to the era 
when |eln \kr^\\ ~ 1. 

We will return to this question in ^ where we will argue that a reasonable 
prescription exists for handling time-dependent logarithms which does not impair our 
ability to extract a predictive initial condition even if |eln|fcr,,|| ~ 1 at the time of 
interest. We can already note one obvious strategy. If only logarithmic divergences are 
present, then better control over the power series may be achieved by "resumming" — 
that is, accounting for the contribution of — all terms at the same order as the leading 
logarithms, which are terms of the form (eln |fcr,,|)"' for all n ^ [53]. For this 
purpose the principal tool is the renormalization group equation. Such an analysis 
has been carried out by Burgess, Holman, Leblond & Shandera [76], whose method 
will be outlined in ^ Unfortunately, this method leaves open the question of whether 
all divergences in the |/cr*| — )■ limit are comparatively tame logarithms, or if more 
aggressive behaviour can occur. 

This question was taken up by Weinberg in 2005 [5l], and later by Chaicherdsakul 
|56j . Weinberg was able to prove that in many inflationary models, but not all, the worst 
divergences would be logarithmic. In a later publication, this theorem was extended 
to fields of higher spin [55]. What of the possibility of faster divergences, which are 
apparently allowed by the analysis of Sasaki et al. [20]? These would grow like powers 
of the scale factor, a(t), corresponding to powers rather than logarithms of \kT^\. In 
many cases such aggressive growth would induce a failure of predictivity long before 
the end of inflation. Also, whereas the coefficient which accompanies a logarithmic 
divergence is physically meaningful, and can be used as an input to renormalization 
group calculations, power-law divergences are by contrast mostly meaningless. The 
late-time physics whose presence they signal must be found elsewhere: it cannot usually 
be distilled from the properties of the divergences themselves. Although Weinberg's 
theorem identifies a class of models containing interactions which may give rise to such 
fast divergences, it does not appear that any previously proposed infiationary model 
makes essential use of such interactions. For this reason they have not yet received 
much attention. 

2.2. Box- cutoff logarithms 

Until non-linear questions became pressing, infiationary perturbation theory was 
dominated by the traditional Lifshitz approach [TTJ [78] , in which one separates each field 
into a background (j) and perturbation 5(j). At second order and above such calculations 
become more arduous, although they have now been carried to an impressive degree of 
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refinement [ISl EHl E] • 

The separate universe principle is an alternative to the Lifshitz approach. According 
to this principle, a volume of spacetime containing a background field and perturbation 
S4> behaves (on scales sufficiently large that gradients may be neglected) just like an 
unperturbed universe containing the homogeneous field value + 6(f). If we solve for an 
arbitrary quantity U with initial conditions set by the value of a background field at 
time T^ in the unperturbed universe, we can determine the values taken by this quantity 
on superhorizon scales in the perturbed universe using the trivial identity 

<5f/(r) = f/(r,0, + (50,)-f/(r,0,), (3) 

where r is the time at which we wish to evaluate 6U. This might typically be at the 
end of inflation or later. Whatever time we choose, Eq. ([3]) allows correlations of 6U to 
be computed provided we know the correlation functions of S(j) at time r*. For example, 
the first contribution to the two-point function oi 6U is {5U{t)5U{t)) = U'^^^{t){6(J)6(J))^, 
where U^^^ denotes the partial derivative of U with respect to 0*. 

Higher-order terms present in Eq. ([3]) imply that other contributions must exist. 
Let us restore spatial arguments for clarity. Even if 6(f)^ has Gaussian statistics at time 
r=K, there will be contributions of the form 

((5t/(r,xi)<5f/(r,X2)) D ^?72^.^,(50(xO50(xi)<5(/.(x2)5(/.(x2)), 

+ gf^,,^.f^,,^.^*0*('50(Xl)(50(x2)50(x2)(50(x2) + (Xi ^ X2)),, (4) 

where (xi f-)- X2) denotes the preceding term with xi and X2 exchanged and the 
symbol "3" means that the correlation function contains this contribution among 
others. Indeed, these four-point correlations will be accompanied by six-, eight- and 
higher 2n-point correlation functions for all n. If 6(f)* has non-gaussian statistics then 
even more contributions will be generated. Zaballa, Rodriguez & Lyth introduced 
a set of diagrammatic ( "Feynman-like" ) rules designed to keep track of these terms 
[82] . Generalized rules were discussed by Byrnes, Koyama, Sasaki & Wands [51] and in 
Ref . [52] , but practical calculations are rarely of sufficient complexity to require them. 
These contributions depend on the correlations among 50^, when |x2 — xi| becomes 
large, which are identified most simply in Fourier space. We find that ((50^(xi)(50^(x2))* 
receives contributions of the form 

^^(xO50^(x2)). D J (0,e-(— ) J ^P(|k - q|)P(,), (5) 

in which P{q) is the power spectrum of 50*, defined by 

(50(ki),50(k2)). = {2n)''6{k, + k2)P{h). (6) 

In the limit |xi — X2I — > 00, significant contributions to the k-integral in Eq. (|5]) arise 
only from the region k ^ |xi — X2|^^. The behaviour of the integrand in the small- 
k limit depends crucially on P{q). For a scale-invariant power spectrum generated 
during inflation P{q) ~ H^/2q^, which gives logarithmic singularities in the q-integral. 
Boubekeur & Lyth [IH], and later Lyth & Rodriguez [SU], remarked that because the 
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logarithmic behaviour is softened in the hmit q ^ k, the log-divergent part only receives 
contributions for q ^ k. Therefore, to a reasonable approximation, 

(502(xO502(x2)). ^ l^—j j — (lnA:L)e''^-(---), (7) 

where L^ |xi — X2I is an infrared cutoff. This expression diverges in the limit L — )■ 00. 

It is clear at once that this divergence is not a physical prediction. It depends on 
what is assumed about correlations among the perturbations in the infinite volume limit, 
where |x2 — xi| — )■ 00. Even within the framework we are using, it is easy to see that 
if P{q) has constant red tilt then this logarithmic divergence is exchanged for a power 
law. If P{q) has constant blue tilt the integral is convergent. An exactly logarithmic 
divergence occurs only for a scale invariant P{q)- For example, if 50 has even a very 
small mass this will give rise to a finite correlation length, beyond which correlations 
decay exponentially. The divergence in Eq. ([7]) would then be absent. We will see later 
that there are good physical reasons to believe this is what should happen in practice, 
at least in certain theories, and some evidence from concrete calculations that it does. 

In any case we should also recognize that the framework we have been discussing 
is inadequate for the description of correlations on very large scales. Why is this? To 
calculate the power spectrum from a model of inflation we must typically assume that 
the background field is spatially homogeneous and depends only on time, whereas 
the spatially dependent fluctuation 50 satisfies |50| ^ |0| everywhere. Although this 
is reasonable on scales not too much larger than the de Sitter horizon, it need not 
happen that every field admits such a decomposition as we pass to the infinite volume 
limit. If this is the case we should set the scale L to be conservatively smaller than 
the largest scale for which a spatially homogeneous background can be found. The 
calculation makes sense within this "box." If we are forced to discuss correlations on 
larger scales, they will have to be determined by patching together a mosaic of boxes in 
which the homogeneous background field may take different values. This point of view 
was advocated by Boubekeur & Lyth [19] , and later refined by Byrnes, Koyama, Sasaki 
& Wands [51], Lyth [83l [52] and other authors [53l [H Eg |86] . 

This can not be the whole story, because Eq. ([71) shows that the mosaicking 
procedure leaves behind terms involving In kL. These logarithms depend on the 
arbitrary scale L. Assuming the physics of the scalar field zero-mode to be quasi- 
classical on very large scales, the missing element is an accurate map of the average 
scalar field value within each box of the mosaic. This map naturally depends on the 
scale L, but the L dependence of physical quantities cancels when we study correlations 
within the mosaic as a whole. These issues will be discussed in more detail in ^ We 
see again that the appearance of infrared divergences is connected with the existence of 
non-trivial background configurations. In the present case the map could consist of an 
inventory of boxes, pairing each box with the average field value within. Alternatively, 
in a theory such as inflation where predictions are essentially statistical, it could simply 
consist of a probability distribution for the average field value, taken over the ensemble 
of boxes. 
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Figure 1. One- loop corrections to tire power spectrum of an active scalar field. In 
the left-hand diagram a loop of spectator scalar fields, represented by dashed lines, 
corrects the two-point function of an active field, represented by solid lines. This loop 
was computed by Weinberg |54( I55j . An error in the numerical coefficient was corrected 
by Adshead, Easther & Lim [58|. In the right-hand diagram, a self-loop corrects the 
same two-point function of an active scalar, first computed in Ref. |57j and again in 
Ref. |58j . In Minkowski space this diagram would factorize, leaving a scale- free integral 
over the loop momentum q. The non-trivial time dependence of de Sitter endows the 
loop with a scale, of the form of Eq. ([5]). This diagram is the leading correction when 
self-loops are included. It would be accompanied by self-loops of the same form as the 
left-hand diagram, which are suppressed by powers of the slow-roll parameter e. For 
spectator fields there is no contribution from the right-hand diagram, so the left-hand 
loop is the leading term. 



2.3. Logarithms from new physics 

In the previous two sections we have emphasized the role of the background field 
configuration in generating infrared divergences: When determining correlations among 
small fiuctuations we take the background to be homogeneous and approximately time 
independent. In practice it may be neither, leading to the emergence of compensating 
logarithms. 

There is another source of logarithmic corrections which is unconnected with the 
background field configuration. Suppose we define some quantum field theory at a scale 
fi. In using this field theory to compute correlations in vacuum — where there is no 
background at all — we are familiar with the appearance of "large logarithms" of the 
form In E/fi, where E is an energy scale characteristic of the correlation in question. 
Corrections of this sort occur in any quantum field theory, and the field theories we use 
to compute infiationary correlations are no exception. Such logarithms were studied 
by Weinberg [51} 155] . who computed the correction induced by loops of N different 
spectator fields in a model where infiation is driven by the vacuum energy associated 
with a single scalar field. 

[See Fig. [H In contrast to Minkowski space, each external leg is evaluated at the 
same late time r^, and therefore these diagrams should be read from the middle to the 
outside. The interior of each diagram functions rather like an instantonqj Breaking the 
left-hand diagram down the middle, four quanta of the spectator field nucleate below 
the de Sitter horizon, and propagate until one pair annihilates at time r' and a second 

f The conventional prescription for calculating correlation functions in an interacting vacuum can be 
related to the Hartle-Hawking state j87j . 
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pair at time t". The product of each annihilation is a quantum of the active field, 

which propagates freely until the surface t = r^ (represented by the outgoing arrows). 

The shared history of these two particles generates a correlation. In an alternative 

and equally acceptable interpretation, a pair of active quanta nucleate. One particle 

propagates directly to the surface t = t^, whereas the second spontaneously fluctuates 

into a pair of spectator quanta at time r', and at a later time t" coalesce to form another 

quantum of the active field before propagating to the final surface. The right-hand 

diagram can be interpreted similarly, with the self-loop either providing a correction to 

the nucleation event, or dressing the propagation of a single quantum on its route to 

the surface t = t^.] 

As in Minkowski space, these logarithms arise from ultraviolet divergent momentum 

integrals which we can study by cutting off the integral for g > A. For a correlation 

dominated by wavenumbers of order k, the typical integrals with which we are confronted 

take a form similar to 

d^g k (k-qY P , A ^ . . , 

^ V H7 _ . In - + finite, (8) 



l^l<yv (27r)3 g^ |q + k|2 Gtt^ k 

where q is the comoving three-momentum which circulates in the loop, and "finite" 
denotes terms which do not diverge in the limit A — ;■ oo|j It is the appearance of 
k in the combination |q + k| which gives the logarithm a non-trivial k dependence, 
and by the usual arguments we can be assured that the coefficient of the logarithm is 
independent of the precise ultraviolet regulator we choose. This does not make Eq. ([8]) 
entirely unambiguous because its interpretation depends on the meaning we assign to A, 
as we shall discuss shortly. Whatever its meaning, we see the emergence of the structure 
E/fi, with E = k characterizing the scale of the correlation, and A = /x representing the 
point at which we define the theory. 

These logarithms are a different species to the infrared logarithms of § §2.1f[2!2| 
which diverged for fixed k when the calculation was taken to occur in a box of infinite 
spatial or temporal extent |§| Eq. ([8]) is finite in this limit, although we will see below 
that quantum loops can also generate infrared divergent terms. It is true that Eq. ([8]) 
has bad behaviour when fc — )■ 0, but a divergence in this limit occurs even at tree- 
levelljlJ We nevertheless include these ultraviolet logarithms in this discussion of infrared 
effects because they alert us to the fact that large-scale (that is, infrared) structure in 

I Other infinite terms may be present, which diverge like powers in the hmit A — > oo, rather than 
logarithms. We are tacitly assuming that these divergences can be absorbed by renormalization-group 
running of the masses and couplings. 

§ Following Weinberg's computation of loop corrections of the form ([8|), and others, Ref. [53] combined 
the effect of secular and box-cutoff logarithms with that of an ultraviolet logarithm at the same order 
of perturbation theory. In this work, the scale A which accompanies the ultraviolet logarithm was 
erroneously identified with the scale L which accompanies the box-cutoff logarithm. This error was 
corrected from v3 onwards of the arXiv version of Ref. |53]. I would like to thank E. Dimastrogiovanni 
for drawing this to my attention. 

II This divergence arises from our convention of normalizing the field modes as Minkowski space 
oscillators deep inside the horizon, but it is not clear any physical effect exists because the divergence 
will be cut off in any model where inflation began at a finite point in the past. Alternatively, Parker 
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the background fields depends on physics at large values of the Hubble rate, H (that 
is, in the ultraviolet). This feature of gravity, which swaps ultraviolet and infrared 
physics, is known from previous applications of the holographic renormalization group 
to domain wall spacetimes which are asymptotically anti de Sitter [951 ESI [97] or de 
Sitter [981 EHl llOOj - We will return to this question when we revisit the mosaicking 
prescription outlined in §2.2[ because it is clear that the large-scale structure of the 
mosaic will depend on what is assumed about the matter field theory in the ultraviolet. 

We must still decide what meaning should be attached to A. In calculating loop 
corrections, Weinberg [5H [55] made use of dimensional regularization. This entails an 
analytic continuation to 3 + e spatial dimensions, after which a limit must be extracted 
by dropping singular terms as e — ?■ jlOlj . Later, Chaicherdsakul [56] and Adshead, 
Easther & Lim [581 1102] employed the same technique. Because A does not appear 
explicitly in this method, being hidden in subtractions associated with the pole at e = 0, 
these papers quoted the loop correction as a multiple of In k but left the scale A implicit. 
Eq. ([8]) used an alternative procedure, imposing a sharp limit on the comoving momenta 
which contribute to each integral. This method was employed in Ref. [57], and later 
by Dimastrogiovanni & Bartolo [59], and gives results in agreement with dimensional 
regularization. 

In either approach, A must be interpreted as a comoving scale, van der Meulen & 
Smit remarked that because comoving scales are not physical, logarithms of the form 
In A/A; were not easy to interpret |103] . It had previously been observed in Ref. [57] 
that if one takes A to be associated with a fixed physical scale Ap, then at an arbitrary 
time the corresponding comoving scale is A = Apa where a is the scale factor. But 
at which time should the cutoff be evaluated? Senatore & Zaldarriaga argued that an 
appropriate choice would be the time of horizon crossing [61], for which a = k/H^ if 
Hk is the value of the Hubble parameter as the wavenumber k crosses the horizon, at 
time T = Tfc. If so we could conclude that Eq. (I8|), and other similar integrals, should 

and collaborators have claimed that a different procedure is required when extracting observables. In 
a free field theory on Minkowski space, the operator product expansion for some local functional 0{x) 
can be defined via normal ordering, 

0{x)0(x') ^ F{x ~ x') -.O^ix): + G(x - x') + ■ ■ ■ , (9) 

where F and G may be singular as x' — > x, and '• • •' denotes terms which are finite in this limit. 
The difference between 0{x)0{x') and the composite operator : 0'^{x) : (which has a finite vacuum 
expectation value) accounts for the singularities of the Wightman functions at coincidence. For an 
interacting field theory on curved spacetime the normal-ordered product :0'^{x) : must be replaced by 
a more general definition of a renormalized composite operator, OJ^{x). Parker et al. propose that such 
operators be defined by an adiabatic subtraction procedure rather than normal ordering. (In principle, 
the curved space operator product expansion of Hollands & Wald [SSI HH ISHl IHIl HI] could be applied 
instead.) In their prescription, contributions to the classical power spectrum should be computed at 
horizon exit using the Fourier modes of C^(a;) rather than 0(x)C'(x'), as is conventionally done. The 
adiabatic subtraction procedure which makes the renormalized composite operator O^ well-defined in 
the ultraviolet also entails subtraction of the low-fc modes which give infrared divergences, at least in 
the massless case j931 194) . At the time of writing, most authors prefer to define the power spectrum 
using the Fourier modes of 0{x)0{x'). 
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be replaced by 

d3, Mk-q)^_ ^^ In A£ + finite. (10) 



(27r)3g3|q + k|2 Qn^ Hk 

A similar conclusion had been reached earlier, by a different method, in Ref. 
This would give a result in precise analogy with the logarithm InE/^ encountered 
in scattering calculations, but seems incompatible with the use of dimensional 
regularization. However, Senatore & Zaldarriaga went on to argue that the free-field 
propagators which appear in the Feynman rules should also be calculated by analytic 
continuation of the Mukhanov equation to 3 + e spatial dimensions. This was not 
attempted in Refs. [511 ESI EEJ |58]. Analytically continuing in this way generates extra 
terms which enable dimensional regularization to reproduce ( TTOl) . At the present time 
it seems unclear whether the horizon-crossing cutoff can be obtained from a more 
fundamental principle, or whether one must simply adopt it as a prescription giving 
reasonable results consistent with the approximate de Sitter symmetryl5ll It is also 
unclear how this form of dimensional regularization should be extended to momentum 
integrals arising from field redefinitions, which are to be discussed below. At one loop 
these can be dealt with using a fixed momentum cutoff, so this ambiguity is not yet 
pressing. 

Corrections to the power spectrum. For convenience, some results on loop corrections 
will be collected in this section. 

Weinberg gave his calculation in the comoving gauge, where the role of active scalar 
field was taken by the comoving curvature perturbation, (. Adshead, Easther & Lim 
worked in the uniform curvature gauge, where the action describing a 3-point contact 
interaction between an active field and a collection of Ng spectator fields s" can be 
written [11] 






O ySaO Saj ^Sq-Sq 



^S^s'X - d"'5^'s'^d's^ ) , (11) 



where d^x = d^xdr, summation of repeated a indices is implied, a prime ' denotes 
a derivative with respect to conformal time and an overdot denotes a derivative with 
respect to cosmic time. In this equation, following the notation of Maldacena [TO] . 
6L/6(f) represents the first-order equation of motion for (p. This term can be removed by 

% Senatore & Zaldarriaga observed that the choice a = k/Hk in the comoving cutoff was equivalent to 
including only quanta beneath the cutoff at the earliest of the times r' and r" in the left-hand diagram 
of Fig. [TJ In a more complicated diagram, one would choose the earliest such time. Essentially the 
same argument was given in Ref. |60j , which discussed the influence of quanta which redshift under the 
horizon at times subsequent to r'. 
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a field redefinitioiiCJ 



^'^ " "" ~ 4M^ \^ ' (s^-^^s^) - -s^s^j . (12) 

The rules for evaluating correlation functions such as (crcr) are discussed in the review 
paper by Koyama elsewhere in this volume [70]. One finds the power spectrum of a to 
be 

where now the time of observation, r*, is chosen to be almost immediately after the time 
of horizon crossing, r* ~ Tfc. 

To obtain the two-point correlation function of 6(f), the effect of the field redefinition 
must be reversed. One finds 

{Scj>{k^)Scj>ik,)), D ((7(ki)a(k2)), 

where (ki ^ k2) denotes the previous term with ki and k2 exchanged, and {AB)i^ is the 
convolution of A and B with argument k. In evaluating this expression we encounter 
three-point correlations of the form (aSaSa) and four-point correlations of the form 
{saSaSpSjs), but the four-point terms give no contribution to ln(Ap//J*). The three- 
point correlations can be evaluated using the general formula for a three-point function 
of scalar fields given in Ref . |11] , yielding 



(15) 



1 IT4 1 

(^(pi)sa(P2)s/3(P3))* = (27r)^5(pi + P2 + P3)jYr^jturir^'' 



-s^ + pM-pI-pD) 



where pt = Pi + P2 + Ps- After some calculation, one finds that the second term in 
brackets (■ ■ ■) does not contribute, and the power spectrum of can be written 

From this calculation we see that field redefinitions reshuffle the coefflcient of ln(Ap/if*). 
Since gauge transformations are a form of field redefinition, we must expect the 
numerical coefficient —3/40 in Eq. ( TT6l) to shift under a change of gauge. In particular, 

"*" This procedure is notationally misleading, because 5L/ 54> is zero by construction (at least to leading 
order) when evaluated on a propagator and therefore gives no contribution whether we subtract it or 
not. The 5L/S(j) term is a proxy for boundary terms, which following the notation of Maldacena [TO] 
have not been written explicitly, which do not vanish on solutions to the equations of motion |11| . It 
was observed in Ref. |104) that the boundary terms and SL/Scj) terms are both removed by the field 
definition, allowing this notational trick to make sense. 
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Figure 2. Graviton loop corrections to tlie power spectrum of a scalar field, calculated 
by Dimastrogiovanni and Bartolo [59]. Unlike the case of scalar loops, the right-hand 
diagram is not slow-roll suppressed compared to the left-hand diagram. As before, the 
interior of these diagrams can be considered as a sort of instanton for the nucleation 
of gravitional (wavy lines) and scalar quanta (straight lines), which propagate to the 
time of observation r* on the external legs of the diagram. 



this will happen if we change to comoving gauge and obtain the correlation function 
of the curvature perturbation, C,. Therefore Eq. ( !T6|) cannot be compared directly to 
Weinberg's result in Ref. [5l], but is an equally good indicator of the magnitude of the 
loop correctionO 

Some time after Weinberg's ^-gauge calculation of this spectator loop [M], the self- 
loop correction to the power spectrum of an active scalar field was obtained [57] U This 
corresponds to the right-hand diagram of Fig. [H and was found to give a correction to the 
power spectrum of the active field immediately after horizon crossing which amounted 
to 

When more than once active species is present one should also consider loops similar 
to the left-hand diagram of Fig. [H with active fields circulating in the interior of the 
diagram. In principle these allow the species of active scalars to fluctuate into each 
other, but have not yet been calculated. 

Corrections from graviton loops. Dimastrogiovanni & Bartolo obtained corrections from 
graviton loops [59J, shown in Fig. [2J Before discussing these, consider the self-loop 
correction for an active scalar field. Wick contraction of an external field operator 
with an operator in the interior of the loop endows the right-hand diagram of Fig. [1] 

* A calculation of this loop in the uniform curvature gauge was given in Ref. 1.58^ , but the contribution 
from the field redefinition was not included. In Weinberg's calculation, given in the comoving gauge, 
a similar redefinition was dropped, which changes the answer. (A numerical error in Weinberg's 
calculation was corrected by Adshead, Easther & Lim [58], but these authors also did not include 
the contribution from the field redefinition.) I would like to thank P. Adshead for communicating the 
outcome of his calculations concerning this issue. 

jl This calculation depends on a set of Feynman rules for the active field, which were obtained in the 
same reference. Unfortunately a sign error was present in one of these rules, which led to an incorrect 
numerical coefficient in the final answer. This sign error was detected and corrected in Ref. 58 , but 
owing to some typographical errors the final numerical coefficient determined by these authors was 
again given incorrectly. The correct coefficient was first given in v3 of Ref. [57] . 
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with a dependence on the external wavenumber. This happens because operators hke 
9~^ apphed to the interior of the loop generate integrals similar to those in Eq. (j8]), 
containing |k + q| in the denominator. In the right-hand diagram of Fig. |2]the graviton 
pair interior to the loop must contract among themselves, and likewise the external 
fields. Therefore no dependence on the external wavenumber is generated, and this 
diagram factorizes into an overall renormalization of the tree-level two-point function, 

where L is an infrared cutoff of the type which appears in box-cutoff logarithms, as 
discussed in §2.21 For the left-hand diagram, one fjud affl 

where "• • •" denotes terms which vanish at late times together with any subleading 
contributions. These examples are interesting because they show that quantum loops 
can generate "box-cutoff" effects in addition to "new physics" logarithms, which was 
suggested but not demonstrated in Ref. |106j . It is not yet clear how to interpret the 
association of a In kL term with a loop containing spin-2 gravitational quanta, which 
lead to traceless deformations of the metric. However, these terms all cancel in the total 
one-loop correction, as we now explain. 

Together, Eqs. (iT6l) . (fTTj) and (iT8l) - (TT9l) give the largest part of the loop correction 
for each active scalar field in a theory containing N^ spectator fieldsljj The terms 
which are neglected by this combination are higher-order in €=„ but do not include 
compensating factors of Ng. First, the /c-dependent logarithms in Eqs. (TT7|) and (TTOll 
cancel between themselves, to leave a logarithm depending only on the ultraviolet and 
infrared regulators. Second, this logarithm cancels with that in Eq. ( ITSl) : only the 
contribution of Eq. flT6l) is left. Therefore, to this accuracy, the total loop correction is 
entirely free of the infrared regulator, L, 



2A;3 V 407r2 M^ H, 



Pr\k) = ^^ ( 1- ^7d^T|lnTf + ■■■ ) • (20) 



Remarkably, this receives contributions from the isocurvature fields alone. It is not 
clear whether this cancellation is accidental, or enforced by some deeper principle. For 
example, in the case of a single minimally coupled scalar field in de Sitter space (obtained 
by passage to the limit e — )■ 0), this absence of one- loop corrections becomes exact and 
is consistent with de Sitter invariance. If the cancellation is accidental there seems no 

ff Dimastrogiovanni & Bartolo quoted their answer in a different form, but Eq. (J19p gives the only 
relevant terms in the limit |fcT*| — > 0. Dimastrogiovanni & Bartolo found a slightly different coefficient 
for each logarithm. As this paper was being completed, a recalculation of the graviton loop appeared 
in Ref. |105| . which agrees with the result given here. (P. Adshead & E. Dimastrogiovanni, personal 
communications . ) 

t Adshead, Easther & Lim argued that for active fields only se/f-loops of the leading-order interaction 
give physical logarithmic corrections ^^ : fo^ the reasons discussed above Eq. P^ . Loop corrections to 
the power spectrum of one active species from other active species are suppressed by slow-roll factors. 
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k-q 




Figure 3. Loop in pure (j>^ theory. A factor of the couphng g is present at each vertex. 



reason why L-dependent logarithms should not appear at subleading order in e, or at 
two loops and higher. It would be of considerable interest to determine whether this is 
the case. 

Burgess et al. [76] emphasized that, if infrared-sensitive contributions occur, one 
should obtain a time independent answer when working with physical ultraviolet and 
infrared cutoffs, Ap and Lp. In their analysis, this was achieved because the ultraviolet 
and infrared terms combined to give InApLp. If one chooses the infrared cutoff to 
be comoving — as we are doing in this discussion — then de Sitter invariance is broken 
and one cannot avoid a time dependent logarithm of the form rejected by Senatore & 
Zaldarriaga as incompatible with eternal inflation. In the present framework, whether a 
similar problem exists depends on the time-dependent distribution of background field 
values which removes L from Eq. (1201) and similar equations. 



(f)^ loop. A somewhat less physical example is the loop in pure 0^ theory, which will play 
an interesting role in the discussion of time dependence in ^ Taking the interaction to 
be 

S,^ jd'xa'hct>\ (21) 

one arrives at the loop of Fig. [31 van der Meulen and Smit calculated this loop, 
accounting for both ultraviolet and infrared effects |103j . It was later recalculated by 
various groups, since it provides a useful test case [HU |60l [76] . One finds 

(j)'^ loop. A loop in 0^ theory was calculated by Tsamis & Woodard |107|ll08j and Petri 
[1U9] . and later by Burgess, Holman, Leblond & Shandera [70]. The interaction is 

54 D I d^x a^ j^5(t)\ (23) 

and the diagram is the same as the right-hand loop of Fig. [T] although with a different 
vertex. The result is 

P*ik) = ^Al + ^,\nkL\n\krA + ..y (24) 
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Other examples of loops have been considered in the literature. In Ref. [50] a set of loop 
corrections to the two- and three-point functions of an active scalar field were calculated 
in a theory where the scalars coupled to the field strength of an Abelian gauge field. 
Although their structure is rather complicated, these loops gave results which are not 
qualitatively different to those we have already studied. 

3. Time evolution 

As in any quantum field theory, the diagrams in Figs. [TH3] are calculated by identifying 
quasi-free propagating modes which travel along the arcs of the diagram and interact at 
vertices. In each calculation discussed above, the quasi-free modes behave like massless 
scalar fields in eternal de Sitter space which have two-point correlations of the form 

(50(ki,ri)50(k2,r2)) = (27r)3<5(ki + k2)^ 



2kf 

^ r (1 - iA;iri)(l + iA;2r2)eH'=i-i-'=2-2) if n < r^ 

^ \ (l + iA;iri)(l-iA;2r2)e^(^2r2-fciri) otherwise' ^ ^ 

where Hi = H{Ti) and ki = |kj|. The magnitudes are equal on-shell, where conservation 
of momentum is satisfied, so fci = ^2- The distinction has been maintained in Eq. fl25l) 
because in practical calculations one must often calculate off-shell where /ci 7^ ^2? only 
taking the on-shell hmit at the end of the calculation [60 | [T5 | [T02| HI] : the same is true for 
higher n-point functions where one must keep the ki all distinct. Eq. (125 p is built out of 
elementary wavefunctions proportional to (1 — iA;r)e''^'^, which asymptotes to unity in the 
infinite future of de Sitter space where r — )■ 0~. Therefore, once correlations have been 
generated using Eq. ( l25l) they persist into the indefinite future. Had we instead taken 
the quasi-free modes to have a mass of order m, we would have obtained elementary 
wavefunctions of the form {—kTy^'^a^^(T)Hu {—kr), where Hi ' are Hankel functions 
of the first and second kind, respectively, and the order of these functions, u, is related 
to the mass m by 

For m ^ H this gives the Hankel functions an imaginary order. In the limit r — > the 
massive de Sitter wavefunction would have an asymptotic expansion proportional to 



~i->=rr"-[-^^^ ' 'r(i + .) -)^--- P^' 

where "• • •" denotes terms which converge to zero more quickly than those which have 
been written. Clearly, correlations are suppressed at late times whenever u < 3/2 
and vanish entirely as r — t- 0^. Having obtained expressions of this form, Riotto & 
Sloth |110j observed that Eq. fl27|) expresses nothing more profound than that a field of 
mass m has a proper correlation length of order 1/m, and in the infinite future of de 
Sitter space the cosmological expansion has carried any two spatially distinct points far 
beyond a proper separation of order 1/m. For this reason, correlations between spatially 
separated points must decay. 
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Weinberg remarked that this mismatch between Eq. f l25p and the decay of 
correlations in the real universe compels us to account for mass terms non-perturbatively 
[55j. In principle we could deal with a term in the Lagrangian of the form m^ScjP' by 
including its contribution in the elementary wavefunctions of the quasi-free modes, or 
by counting it among the interactions. If we include it among the interactions, then we 
compute correlations using Eq. fl25|) and are obliged to make the implicit assumption 
that they do not subsequently decay. This may yield reasonable results, supposing the 
correlations evolve only slowly, provided we do not ask for their properties too long after 
horizon crossing. Nevertheless, where the correlations truly decay this approximation 
must lead to trouble at late times. Since we expect most fields to acquire a non-trivial 
potential except where protected by an exact symmetry, we also expect the fluctuations 
in these fields to develop small masses. As was explained in §2.11 it is our failure to 
account for this time evolution in the background fields, and the subsequent decay of 
correlations, which is responsible for the appearance of the secular logarithms In |/cr^|. 

What if a field has a potential containing no mass term, but other higher-order 
interactions? This situation was analysed by Senatore & Zaldarriaga [61], who argued 
that on any potential of this kind the background field would still roll down and evolve. 
For this reason, we would expect a mass term to be generated for the fluctuations after 
resumming a sufficient number of insertions in the propagator. This would subsequently 
cause correlations to decay. Therefore this situation is qualitatively the same as the case 
of a mass term in the potential. We will see an explicit example in §3.31 

3.1. Weinberg's theorem and related results 

How fast can we expect correlations to evolve? In §2.11 we discussed a theorem due 
to Weinberg ^54j which guarantees that only logarithms of |/i;r^,| appear in inflationary 
correlation functions, and not power laws of the form |/cr,= |~" for n ^ 1 which in principle 
are allowed [20]. (Contributions proportional to positive powers of \kT^\ are generically 
present, but contribute nothing at late times.) The precise statement of this theorem 
applies to scalar and tensor excitations. In three dimensions, we have: 

Theorem. In three space dimensions, power law divergences are absent at late times in 
correlations among scalar and tensor quanta, provided all interactions in the Lagrangian, 
when written in cosmic time, fall into one of two classes: 

(a) Interactions containing strictly less than one factor of a{t). 

(b) Interactions which may grow as fast as a{t) (but no faster), and which contain only 
fields rather than time derivatives of fields. 

Weinberg proved this theorem by studying commutators among the fields and their 
derivatives at late times, and then counting how these commutators could appear in the 
possible elementary interactions. This method of proof does not make use of the explicit 
slow-roll, massless solution (!25|) . However, Senatore & Zaldarriaga later remarked that 
a refinement of this theorem might be possible [61] , because the asymptotic estimates 
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obtained by Weinberg did assume that the fields were massless. In addition, since 
each elementary interaction was considered in isolation, the possibility of cancellations 
among groups of interactions was not considered. At the time of writing, these possible 
refinements have not been studied. 

In a later publication, Weinberg extended this theorem to Abelian and non-Abelian 
vector fields, and Dirac fermions [55j (see also Ref. ^6]). These fields have rather better 
behaviour than the scalar fields considered in the original theorem, and correlation 
functions among them typically converge at late times. 

3.2. Resummation by the SN formula 

For the remainder of this discussion, we will restrict attention to theories satisfying 
the conditions of Weinberg's theorem. Therefore only secular logarithms will appear as 
T^ — 7- 0, and not the faster power-law divergences. 

Up to this point we have not made our choice of gauge explicit, but we have been 
discussing fiuctuations in scalar fields. This will usually correspond to the uniform 
curvature gauge. The quantity whose correlations are accessible to observation in the 
CMB is often taken to be the comoving curvature perturbation, (. (Indeed, Weinberg's 
theorem and the spectator loop discussed in §2.31 were given directly in this gauge.) 
The statistical properties of ( are determined from a superposition of the active fields 
in the theory. Many prescriptions have been employed to obtain this superposition 
|llll 1112] , but a simple method is to observe that C is a local fiuctuation in the aggregate 
cosmological expansion, ( = 6lna, where a is the scale factor. We can therefore use the 
separate universe formula, Eq. ([3]). One usually defines Ina/ao = A^, where A^ is the 
number of e-folds from a reference time where a = a^. We take the initial slice to be 
a uniform curvature hypersurface, and the final slice to be a surface of uniform energy 
density, after which we arrive at the 6N formula |113l 11141 150] , reviewed elsewhere in 
this volume by Tanaka, Suyama & Yokoyama [115] and Wands |116j 

Among other things, this is a gauge transformation. The 6N formula reproduces those 
terms from the non-linear transformation between the uniform curvature gauge and 
comoving gauge [TU] which do not vanish on superhorizon scales. 

Eq. fl28|) enables the correlation functions of ( to be determined from those of 
6(j)^:, which in turn are computed using the massless propagator ( 12^ . In the foregoing 
discussion we argued that Eq. ( 125!) prevents correlations from decaying after they have 
been generated, and that this corresponds to neglecting evolution of the background 
fields. The 6N formula ameliorates this difficulty. Suppose we wish to compute 
correlations on some scale k. If we take the initial time, r*, to be shortly after the 
time of horizon exit associated with k, then we do not commit a gross error in failing 
to account for their subsequent decay. Moreover, the time dependence summarized by 
A^ accounts for evolution of the background fields as they roll down their potentials. 
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Therefore we expect Eq. (125]) to incorporate the infrared physics whose neglect gave rise 
to secular logarithms. This argument was given in Ref. [53] and is supported by the 
analysis of van der Meulen & Smit |103] , although a formal proof is still lacking. 

The same conclusion can be reached by considering the logarithms themselves. If r* 
is chosen not long after horizon crossing, then the terms In | kT^, \ are of order unity. If the 
leading term containing n powers of In |A;r^,| occurs in the combination ~ (eln |A;r^,|)", 
then, not long after horizon crossing, all terms containing logarithms will be small 
compared to the tree level whenever e ^ 1. Therefore we can neglect all logarithms in 
using Eq. (128]) to compute correlations of (. The effect has been to "resum" their effect 
into the background evolution of the scalar fields, which is described by the accumulating 
number of e-folds, N. For example, in a single-field model of inflation this argument 
reproduces the expected conclusion that ( does not evolve. But the virtue of the SN 
method is that it is not restricted to the single field case: In a multiple field scenario 
where ( has non-trivial evolution, the SN formula will correctly describe the growth and 
decay of fluctuations. 

There may be problems when using Eq. fl28]) to obtain the correlation functions of ( 
if one wishes to include the "quantum" logarithms of §2.31 This procedure was applied 
in Refs. [STJ [531 EH EH] [60]. However, as we have already observed, the 6N formula is a 
gauge transformation, and it was argued in §2.3l that such transformations will typically 
modify the coefficient of In Ap/Mp. Although this possibility has yet to be investigated 
in detail, the large number of ultraviolet terms discarded by the 6N formula (which 
should be compared, for example, with Eq. (A. 8) of Ref. [lO]) would apparently make a 
mismatch likely. The 6N formula retains all relevant infrared physics, so we would not 
expect similar difficulties with secular and box-cutoff logarithms. 

3.3. Resummation by the dynamical renormalization group 

For the purpose of comparison with observation, the 6N method is presumably sufficient 
and allows us to extract initial conditions from typical inflationary models without any 
loss of predictivity. In many models of inflation we wish to obtain the properties of 
correlations only 50 to 60 e-folds after they are synthesized. In a conventional model 
the only source of appreciable time evolution during these 50 to 60 e-folds comes from 
the classical variation of the background. Nevertheless, as a point of principle, one might 
prefer a more satisfactory resolution. Some time after the resummation interpretation 
of 6N had been proposed. Burgess, Holman, Leblond & Shandera remarked that it did 
not account for the possibility of secular time evolution from quantum effects [76], such 
as those of Eqs. f[T9]) . f l22]) and 1^^. These terms typically appear in combination with 
the small parameter {H/Mp)"^, and are therefore small unless many e-folds have elapsed 
since horizon crossing. 

Burgess et al. suggested that these logarithms could be resummed using the method 
of the renormalization group, and applied this technique to the quartic loop of Eq. ( HM . 
We briefly review their method. The starting point is to introduce an arbitrary constant. 
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c, which is unity in the tree-level power spectrum, 

n(c,fc) = ^. (29) 



Our strategy is to resum the large logarithms into c. To do so, rewrite Eq. flMj) in terms 
of an arbitrary intermediate time scale "d, 

Clearly, the combination of overall prefactor and first bracket is of the same form as 
P^{c,k), with c = c('i9). We can therefore replace one by the other. Whatever the 
result, it must be independent of the arbitrary time "(9. Applying the familiar method 
of Gell-Mann and Low [117J, we find 

Solving this differential equation and setting {} = \kT^\, we conclude 

PM = ^\kn\' (1 + 0[eMn |A;r,|]) , (32) 

where S satisfies 

6 = -^^\nkL. (33) 

The error in this procedure is of order e^ln|A;r*|. Therefore Eqs. fl32|) - fl33|) are 
trustworthy even when |eln|fcr*|| ~ 1, giving rise to a resummation of the leading 
logarithms. Subleading logarithms could be included by accounting for terms of order 
A^ or higher in Eq. (I3T]) . 

Eq. (152]) describes decaying correlations. On the basis of the foregoing discussion, we 
expect this behaviour to represent the correct infrared physics in a large class of theories. 
Indeed, comparison with Eq. fl271) shows that the effect of c(|fcr^.|) is to introduce a 
dependence on |fcr*| comparable to the case of a massive scalar field. Making use of 
Eq. (l26l) and recalling that the square of Eq. (1271) yields the power spectrum, we see 
that Eq. (15^ is equivalent to a dynamically generated mass M^s [75] . 

M^s = ^H'jnkL. (34) 

If 1/k and L do not generate an exponential hierarchy and A ~ 0(1), then Mes is of 
order the Hubble rate. We can regard this as a concrete example of the argument of 
Senatore & Zaldarriaga [^I] discussed just before §3.11 

This conclusion does not apply to every theory. Burgess, Holman, Leblond & 
Shandera observed that the same argument applied to pure 0^ theory does not yield a 
dynamically generated mass. The obstruction comes from the power of In |A;r^,| in the 
leading term of Eq. (122]) . Ignoring the ultraviolet logarithm proportional to In Ap/Mp, 
which plays no role here, resumming the secular logarithms would give 



P.ik) = gexp (^^^Mkr.ir) {l + 0[eHn\kr^\ 



(35) 
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Since this does not suppress the power spectrum by a positive power of |A;r,,| at late 
times we cannot interpret it in terms of an effective value of u in Eq. fl27j) . and thus 
an effective mass Meg. Indeed, ( l35l) has rather poor behaviour in the limit \kT^,\ — )■ 0, 
because the exponential diverges there. One should not take this behaviour literally. 
Burgess et al. included logarithms of cubic order which have been neglected here, and 
which would become negative in the limit r* — )■ 0^. Such terms would compete with 
the positive contribtion from (ln|fcr*|)^. Thus, in pure 0^ theory the behaviour of the 
exponential in the far future becomes a delicate question. 

This is not difficult to interpret. In pure 0^ theory the Hamiltonian is unbounded 
below, and there is nothing to prevent fluctuations (and their correlations) growing 
indefinitely. Burgess et al. remarked that this pathological behaviour made studying 
the late-time behaviour of this theory problematic. Following the argument of Senatore 
& Zaldarriaga [61], we may expect that in realistic theories behaviour comparable to ( l32l) 
is more generic. 

4. Spatial evolution 

Less is known about the resolution of what we have called box-cutoff logarithms, of the 
form In kL, which were left behind in §2.21 after tiling the spatial region of interest into 
a mosaic of boxes. 

4.I. Mosaicking prescriptions 

Our discussion in §2.21 was based on an assumption: that the role of physics at scales 
larger than L was to determine a quasi-classical background field configuration in 
each box. A proof of the validity of this assumption is not known. What would 
it entail? Consider an expectation value of local operators Oi, . . . , 0„ evaluated at 
spatial positions xi, . . . , x„ with roughly common separations |xj — Xj| ~ 1/k which are 
small on the characteristic scale of the mosaic, so that kL ^ 1. We wish to evaluate this 
correlation function at time t^,, just after horizon exit of the wavenumber k, which can be 
achieved by averaging against some wavefunctional "^[O] weighting spatial configurations 
of theOi [2n|fTT8]. 

(Oi(xi) ■ --OnM)* = f[dO] Oi(xi) ■ ■ ■0„(x„) |^*[0]|2, (36) 

where a subscript '*' denotes evaluation at t^. Eq. f l36|) leads to the in-in or Schwinger 
path integral once each wavefunctional has been written as an integral over all histories 
0{t, x) with appropriate boundary conditions at r*. Cut out a box of comoving size ~ L, 
approximately centred on the Xj. The path integral divides into separate integrals 
over the spatial field configuration interior to the box, written 0~ , and the exterior 
configuration, O"*". These configurations merge continuously on surfaces of the box, 
denoted dL. In these terms, Eq. (1361) can be rewritten 

(Oi(xO • ■■O^i^n)). = J[dO-dOg^] Or(xO • • ■0„(x„) |K[0-,0,J|l (37) 
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The interior wavefunctional \l/~ satisfies 

K[0-,0,^] = J[dO-{T)] exp {iS[0-]) y"[dO+(r)] exp (i^[0+]) , (38) 

where S" is a local action functional. In this expression, [dO,J'(r)] represents an integral 
over interior field histories 0~{t, x) which coincide with 0^(x) at t^, and have boundary 
data Oq^ on dL, whereas [dO"'"(r)] represents an otherwise unconstrained integral over 
all exterior field histories 0'^{t, x) which merge sufficiently smoothly with the boundary 
configuration Og^ on the edges of the box. 

In typical applications, the expectation value ( l36l) does not depend on the prior 
history of the universe up to the point of horizon crossing. It is therefore independent of 
almost all boundary data represented by Oqj^, retaining only a possible dependence on 
the boundary configuration near the surface t = t^,. If the path integral is dominated by 
semiclassical field configurations which vary only slowly on the scale L, then to a good 
approximation this remaining dependence will involve only the roughly homogeneous 
boundary value at r*, and not its spatial gradients. Eq. f p8|) suggests that expectation 
values are to be calculated by fixing a homogeneous background field configuration at the 
surface of horizon crossing, and summing over unrestricted field histories for fluctuations 
in the interior. Finally, one averages over the fixed background configurations with an 
unknown probability measure, P, which derives from the integral over exterior field 
configurations. At this point we have arrived at the prescription of §2.21 

This construction makes P a distribution over field configurations at a fixed spatial 
position in an ensemble of realizations, but if ergodicity applies we can equally interpret 
it as a probability distribution over different boxes in the same realization. The above 
argument is merely heuristic, and the existence of such a decoupling limit has not yet 
been demonstratedjjj Nevertheless, Eqs. f p6|) - fl38|) show a strong similarity with path 
integral approaches to the operator product expansion ( "OPE" ) in fiat space quantum 
field theory |119] . A manifestly local and covariant formulation of this expansion in 
curved spacetime has been given by Hollands & Wald [HI [90l [921 [89l [88], which in 
principle could be applied to the study of infrared effects in nearly de Sitter spacetimes. 
However, this has not yet been done. 

A similar conclusion was reached by Allen & Folacci |120] and later Kumar, Leblond 
& Rajaraman |121] . who studied ambiguities in defining the zero-mode for the Bunch- 
Davies propagator associated with a massless, minimally coupled scalar field in de Sitter 
space. Kumar et al. noted that similar ambiguities were removed in Minkowski space 
owing to requirements imposed by the principle of cluster decomposition. Exploiting 
the extra freedom available in de Sitter space, they succeeded in subtracting infrared 
divergences by performing an appropriate redefinition. In a field theory this redefinition 
would naturally become spatially dependent on scales which are large compared with 
those appearing in expectation values. The result is to reproduce the conclusion that 

f If the probability distribution on the boundary does not decouple from spatial derivatives, or the 
field configuration prior to horizon crossing, then the mosaic approach wc have adopted need not apply, 
and the whole question of infrared effects becomes more interesting. 
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infrared divergences derive from a supposition that calculations are to be carried out with 
a definite, homogeneous expectation value for the background field configuration. From 
this point of view, Eq. ( l38l) can be thought of as an analogue of cluster decomposition 
in de Sitter. 

These ideas have a large literature of their own. Kirsten & Garriga jl22] used related 
reasoning to represent the quantized massless scalar field as the product of a Hilbert 
space, associated with the zero-mode, and a Fock space corresponding to excitations 
of finite wavenumber. Related ideas were explored by Moncrief |123] . Higuchi and 
collaborators [IH lESl 1261 [I27] , and later Giddings & Marolf [I28]. Urakawa & Maeda 
investigated the same physics from a different perspective [1291 1130] . 

4-2. The L-dependence of physical predictions 

While discussing logarithms generated by ultraviolet effects, Senatore & Zaldarriaga 
remarked that combinations such as InA/fc in Eq. (E]) were unphysical [61], because 
there is no physical comoving scale which could accompany k. However, this argument 
does not entirely preclude the appearance of k in any logarithm. For example, the 
combination ln|fcr^,| measures by how many e-folds the wavenumber k is outside the 
horizon at time r* and is a sensible physical quantity. Nevertheless, this example shows 
that one must be careful in dealing with the combination In/cL, in which L is also a 
comoving scale. 

Let us first idealize to an inflating volume which is spatially infinite. As in §2.21 we 
divide this volume into boxes of size L, in each of which the field is spatially homogeneous 
with small perturbations. We have seen that correlations computed within these boxes 
depend on the arbitrary scale L, but (assuming only zero-modes need be retained from 
the long wavelength physics) an answer independent of L is obtained if we average the 
correlations over boxes |131j . This requires a knowledge of the homogeneous background 
field within each box of the mosaic. If merely probabilistic answers are acceptable, it 
may only be necessary to use information about the statistical distribution of field values 
over the ensemble of boxes. Cosmological perturbation theory has been discussed from 
this point of view by many authors; see, for example, the recent textbook by Mukhanov 
[132]. 

If we assume that only zero-modes need be retained, as in §4.11 it follows that this 
argument does not depend on a careful discussion of how the total L dependence can 
somehow cancel between the L-dependent box-sized correlations and the distribution of 
field values among boxes. The schematic argument of Eqs. (!36l) -( l38l) shows that in a full 
treatment, where the entire ultraviolet physics and initial conditions are known, we have 
done nothing more than divide the calculation into two parts. When we assemble these 
two halves, the final answer must be independent of L. This happens automatically, 
no matter how remarkably fine-tuned the necessary cancellations appear, and is a 
standard argument in the application of effective field theories [72] . The scale L therefore 
disappears from the answer and we need not worry about its physical significance. In 
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practice, the distribution of field values among boxes of the mosaic will evolve in time. 
After convolving with this distribution, L-dependent terms would be traded for a time 
dependence which could itself be represented by secular logarithms. 

However, it may not be the case that this L-independent answer is the quantity we 
wish to compare with observation. In discussing similar questions almost twenty years 
ago, Salopek & Bond remarked that we only wish to calculate what we can observe 
|133[ 1134] . What can we observe? Only the density fluctuation within our presently 
observable patch of the universe |135l [83| l52l [8il 185] . Observationally it is a matter of 
perfect indifference to us (although, of course, of surpassing theoretical interest) whether 
the distant, unobserved universe contains many regions similar to our own, or whether 
all distant regions are quite dissimilar. Whichever is correct, the answer will not change 
the result of satellite observations made today. 

If we subscribe to programme outlined in ^D then we must suppose that our 
observable region of the universe passed through a phase of inflation which eventually 
came to an end. We assume that the classical background associated with whatever 
scalar fields supplied the necessary vacuum energy became roughly homogeneous within 
a box somewhat larger than our present horizon. This box should be large enough 
to contain many regions of comparable size to the presently observable universe, but 
need not be exponentially larger [831 El]. We wish to compute correlations in the 
density fluctuation which would be recorded by a typical observer whose local patch 
experienced the same gross cosmological history. For this reason it may be incorrect to 
compare the L-independent answer discussed above with experiment, because in this 
answer the correlations specific to our patch are commingled with correlations recorded 
by observers experiencing a gross cosmological history quite different to our own. 

Therefore, let us choose L so that In kL ~ 1, making L a little larger than our 
presently observable universejjj The problem we face in this scenario is to determine 
the correct homogeneous background in which we should carry out our calculation. In 
general the answer depends on the distribution of field values within the ensemble of 
boxes, but there is a simple class of models in which this question is trivial. In any single- 
field model with a unique reheating minimum there is a unique value of the background 
field, = 060; when inflation is of order 60 e- folds from ending. The prior history of 
the large-scale universe is irrelevant [251 ^^. To compute what would be recorded by 
a typical observer in this model, we are entitled to carry out our calculation within 
a box carrying an approximately homogeneous background field with value ~ 06o- 
In this picture, the effective field theories we use to study slow-roll inflation are to 
be understood only as a description of the process by which an inflationary trajectory 
arrives in its final reheating minimum, causing inflation to end. 

f This prescription seems to have been rediscovered several times, beginning with Salopek & Bond 
[1341 I133J . In our present context the choice In fcL ~ 1 was suggested by Boubekeur & Lyth [49J , who 
referred to it as a 'minimal box.' Its properties were later studied in a series of papers by Lyth and 
collaborators [5D1 |HS1 113 IHSj ■ The same prescription, in various forms, was given by Bartolo et al. [Hi] 
and later by Enqvist et al. '85] and Kumar et al. [136] . 
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In the distant future our observable region of the universe will expand to include 
modes which are presently unobservable. Enqvist, Nurmi, Podolsky & Rigopoulos 
remarked that an inflationary theorist, calculating millions of years in the future, would 
be obliged to work in a box of slightly larger size V > L [85j. The argument of Boubekeur 
& Lyth |19] , which was extended to two loops by Bartolo et al. ^84j , shows that physical 
observables compose correctly under the process of averaging L-sized boxes to make an 
L'-sized box. 

4-3. Stochastic inflation 

In a scenario which is more general than single-field inflation there may be many 
terminal vacua in which the universe can reheat, rather than a unique choice. In these 
models there is no preferred background configuration — analogous to the single-field 
configuration ~ (/)go — in which we should carry out our computations. Even where a 
terminal vacuum can be selected in advance, it was shown by Byrnes, Choi & Hall that 
the final fluctuations may depend sensitively on the path by which the fields arrive at 
the minimum [1371 1138] (reviewed in Ref . |139 ] ) . 

To determine what would be recorded by a typical observer we must know at 
least the frequency with which these minima are populated by boxes in which inflation 
terminates. This can be ascertained only if we know the distribution of field values 
within boxes of the mosaic, and (as we have discussed above) this distribution is an 
ultraviolet-dependent quantity. Without knowledge of the complete ultraviolet physics, 
including any relevant initial conditions, we cannot calculate it from first principles. 

Apart from our interest in observational predictions, there may be other reasons to 
enquire about the large-scale disposition of the scalar fields. Senatore & Zaldarriaga [61] 
emphasized that we would like to know the circumstances under which eternal inflation 
can occur |140l 1141] . It is possible this depends on effects arising from loops at high 
order in perturbation theory. For this reason we may wish to go beyond the narrow 
view advanced above, according to which our theories of slow-roll inflation should be 
restricted to a description of how inflation ends. Using a model of slow-roll inflation 
to describe evolution within a large spatial volume over very long times would again 
require us to know the distribution of field values over the entire mosaic of boxes. 

Clearly we cannot hope to calculate this distribution ah initio. Suppose, however, 
that by some means we know the distribution of field values over the mosaic of boxes 
at time tj, at which point we suppose the Hubble scale took the value Hi. If we have an 
effective description of the processes by which correlations are established at energies 
lower than H^, then we can calculate the subsequent evolution of the mosaic. 

Similar techniques are used in many applications where the long-range physics 
likewise cannot be determined from first principles. A simple example is the distribution 
of partons within colliding nuclei |142[ 1131] . Since this distribution depends on the 
unknown details of confinement in QCD it cannot be calculated from first principles. 
However, its evolution with energy can be calculated at energies large enough that 
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the perturbative regime of QCD is a good description. In the inflationary case, it is the 
regime of smaU H/M-p which is accessible to perturbative calculations. A simple method 
of calculating this evolution was suggested by Starobinsky |143] and later applied to the 
calculation of correlation functions in the deep infrared by Starobinsky & Yokoyama 

In Starobinsky's method we suppose that the background field value within a given 
box of the mosaic is known. Our discussion of time evolution in §2.11 and [|3] shows that 
we will obtain reasonable answers only if we account for the slow roll-down of this field. 
In a single-field scenario the evolution of the background field value is given by 

V' if3/2 

+ -l^O{t). (39) 



2,H 2ti 
where V{(j)) is the potential and 9{t) is a Gaussian noise term satisfying 

{e{t)e{t')) = 6{t-t'). (40) 

In a careful treatment, we would find that this noise term depends on what we assume 
about the physics responsible for generating correlations within this box. Eq. fl39l) is a 
Langevin equation, which can be equivalently expressed as an evoution equation for the 
distribution function of field values over the mosaic, /(0), 

df _ 1 difV) H-' d'f 

dt 3H d<f) 87r2 902 ■ ^ > 

Various refinements of this equation are possible, which were discussed by Salopek & 
Bond [T33] . 

It has been suggested that stochastic infiation may represent one possible method of 
determining the distribution of field values over the entire mosaic of boxes [SU |85l 1131 



Although this idea is attractive, there are several drawbacks. First, Eqs. (l39ll -( HT1) 
depend on an initial condition for the distribution f{4>). This initial condition is sensitive 
to the details of ultraviolet physics and is no more calculable from first principles than is 
/(0) itself. This difficulty is usually circumvented by seeking stationary solutions of the 
Fokker-Planck equation, Eq. (j^Tl) . However, it is not clear under which circumstances 
this is a reasonable choice, or whether the result genuinely approximates what would be 
obtained if we were to retain the complete ultraviolet physics. 

Second, the evolution equations (J39ll and (l4T1l depend on what we assume about 
physics at high values of the Hubble parameter, up to if ~ Mp. In particular, these 
equations assume that slow-roll infiation is a good description throughout this regime. 
This is analogous to the supposition that the Standard Model is a good description of 
physics at energies from the electroweak frontier at ~ 1 TeV to up the GUT scale at 
~ 10^^ GeV. This assumption allows renormalization group evolution of the SU(2), 
U(l)y and strong coupling constants to predict a unification near the GUT scale in 
the same way that Eq. (j^Tj) allows us to predict the behaviour of (p ^'^ values of H 
far above those which gave rise to observable density fiuctuations. In either case it is 
perfectly reasonable that unguessed physics intervenes to spoil the conclusion. If so, 
the true evolution of the background fields might follow laws very different to Eqs. (!39ll 
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and f l4ip . Alternatively, at sufficiently high values of H the concept of a homogeneous 
background field itself may lose meaning. In that case the present discussion would be 
invalidated, becoming relevant only at lower energies. 

It was emphasized above that, assuming only zero modes need be retained from 
the quasi-classical physics on large scales, cancellation of L was automatic after 
averaging over the distribution of field values within the mosaic. Since there is no 
recipe to construct this distribution, it must checked individually for each candidate. 
Unfortunately, it is non-trivial to verify that the L-dependence carried by a solution to 
Eq. (14T!) would precisely cancel that L-dependence of correlations computed within a 
box L. At the time of writing, this does not appear to have been verified. A complete 
cancellation would depend on retaining sufficient information about the relevant long- 
range physics in Eqs. (l39l)-( l4T|) . and it is far from obvious that this is the case. However, 
since the same physics which gives rise to the box-cutoff logarithms is used to compute 
the evolution of in Eq. ( 139|) . it is at least plausible that the L dependence of / may 
have the requisite properties. A similar property is exhibited by the parton distribution 
functions which describe the physics internal to colliding nuclei |142j . A ffist step 
in this direction was given by Kuhnel & Schwarz |145j who argued that large-scale 
correlations were strongly suppressed by stochastic effects, leading to infrared finite 
correlation functions. 

4.4- Dynamically generated masses 

In §2.11 and ^ we discussed the role of an evolving background field in suppressing 
correlations at late times, and in §3.31 we reviewed the argument of Burgess et al. 
[76] , according to which the secular logarithms could be resummed into a dynamically 
generated mass. This mass drives any correlations to zero in the far future, as the field 
settles in its minimum and loses memory of its past history. 

Eq. fl39l) and its associated Fokker-Planck equation, Eq. fj4Tl) . account for both 
fluctuations and time evolution in the background field. A fluctuation 9 in Eq. ( l39l) 
can push the field up the potential. At later times the term —V'/3H causes it to roll 
towards the minimum. Therefore we would not expect statistical properties computed 
using the stochastic framework to suffer from the pathological failure to decay associated 
with correlations established using Eq. fl25|) . This expectation is borne out in explicit 
calculations. 

Starobinsky & Yokoyama worked with a model of A0^ inflation, defined by the 
interaction of Eq. (l23l) . and calculated the evolution of correlations in this theory using 
the stochastic method [144j . Working in the massless case for simplicity, they were able 
to use the Fokker-Planck equation fHTl) to derive an evolution equation for the two-point 
function (0^), 
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Starobinsky & Yokoyama observed that the solution of this equation gave a constant, 

<*^> - IS <«' 

at late times, and remarked that this was equivalent to the dynamical generation of a 
mass of order H. We briefly review this conclusion using the analysis of Riotto & Sloth 
PUJJ . who derived a "gap" equation for the (j) propagator, G, 

(d + - \(4>^) + G(x, x')] ] G(x, x') = -^6(x - x') (44) 

where the differential operator D = g°'^'VaVb operates on the variable x, and Va is 
a covariant derivative compatible with the metric Qab- Eq- (Hll) is obtained using 
the method of the 2-particle irreducible action. For details the reader may consult 
the original article by Riotto & Sloth |110j . or the textbook by Calzetta & Hu |146j . 
After translating to Fourier space and substituting Starobinsky & Yokoyama's late-time 
solution for (0^), given by Eq. fj43|) . in Eq. dH]) one finds that the the Green's function 
G{x, x') should be built out of solutions to the homogeneous equation 



,2 47,- 



G + "iHG +\^ + ^H^ G = 0, . (45) 



This is the Klein-Gordon equation for a field of mass 

M^, = ^H\ (46) 



Later, Burgess et al. [76] argued that this mass could be reproduced from the dynamical 
mass generated by the renormalization group, Eq. flMl) . at least in the large- A^ limit of 
an 0(A^)-symmetric theory similar to that studied by Riotto & Sloth |110] and Petri 

[Tug] . 

4-5. Tensor-to-scalar ratio 

If an exact shift symmetry protects the potential of some field, then it experiences no 
classical evolution. An example might be a Goldstone boson or axion, or more generally 
any isocurvature field which can be interpreted as an angle. Once perturbations have 
been generated in such a field, they do not decayqj In this case we would not expect 
resummation to generate a mass, and there is no prediction for the value taken by such 
a field within our presently observable patch, although it may be possible to estimate 

f During inflation, quantum fluctuations will completely randomize the value of such a field- 
Nevertheless, the existence of a shift symmetry means the field value is unobservable and there seem no 
unacceptable consequences associated with this behaviour. If the field is massless but interacting — such 
as at a critical point — then the field value is in principle observable. (This example was suggested by 
C. Burgess and L. Leblond.) For sufficiently strong interactions the field presumably rolls down the 
resulting potential at late times, as described by Senatore & Zaldarriaga [61j . preventing an unphysical 
divergence to arbitrarily large field values. However, it is not known what conditions must be satisfied 
by the interactions for divergences to be absent. 
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a distribution [H21 ES] • Hertzberg, Tegmark & Wilczek |147] argued that such scenarios 
represented a test case of anthropic reasoning. 

What would happen if a microwave background observable depended on the value 
of such a field? We would not be able to predict which measurements should be recorded 
by CMB survey satellites, such as the WMAP and Planck satellite missions. A similar 
argument was given by Sloth, who pointed out that an observable such as the tensor 
fraction, r, may depend on a ratio of two quantities which receive loop corrections scaling 
differently as the box size increases |l6l H?]. The expected value of r will depend on 
the box in which correlations are computed. However, as has been stressed throughout, 
to obtain observables one must compute a conditional expectation value, contingent on 
inflation ending in our own vacuum (almost surely) everywhere within the box. For this 
reason, the observational relevance of infrared corrections in r is not yet clear. 

5. Discussion 

Infrared effects have attracted considerable interest, but at the time of writing it has not 
been possible to extract a definitive prediction of a novel effect from these calculations — 
that is, an effect which is not already implicit in calculations carried out locally around 
background fields set at a point on the scalar potential. 

The best studied examples of infrared effects are the time-dependent secular 
logarithms introduced in §2.11 Over short timescales these logarithms represent 
evolution of correlations, as the background fields slowly roll down their potentials. 
At later times quantum effects can become important, but their role remains somewhat 
unclear. The calculation of Burgess et al, supported by the arguments of Senatore & 
Zaldarriaga, suggests that in stable theories they will suppress correlations at late times. 

Box-cutoff logarithms represent a different problem: the possibility that, after many 
e-folds of inflation, spatial gradients develop in the background scalar fields. If we 
only wish to use our theories of inflation to study observable correlations generated on 
approach to reheating, then this may not be a significant problem. On the other hand, 
if we wish to study the evolution of a large spatial volume undergoing inflation for a 
long period of time, we expose ourselves to the possibility of unknown new physics at 
large values of H/M-p. The gauge invariant arguments of Losic & Unruh appear to 
cast doubt on the idea that perturbation theory remains globally well defined for an 
indefinite period in de Sitter or nearly de Sitter spacetimes. 

Except in special cases, such as single-field inflation, the study of correlations 
generated near exit from inflation is still ultraviolet-sensitive (in the sense of physics 
at large values of H/M-p) and infrared-sensitive (in the sense of physics on very large 
lengthscales) , because we do not know in which reheating minimum inflation is likely 
to terminate. This need not be a problem if observation can accurately determine the 
properties of our local minimum. In that case we would be able to understand all the 
properties of our vacuum, but we would not be able to understand how to embed this 
minimum in the larger landscape of the inflationary potential. However, this problem 
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is not new and has been with us since the early days of the inflationary paradigm. 
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